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Abstract 

Multiwalled carbon nanotubes based nanomachines driven by thermal bath are proposed in this 
Letter. Atomistic simulation shows that intershell translational and rotational motion can be 
activated at sufficiently high temperature, where the van der Waals potential barrier is accessible 
by the thermal fluctuation. The thermal fluctuation is further identified to be consistent with the 
equipartition theorem. This model can be used to construct one dimensional devices driven by 
thermal bath or nano-channels for molecule transport. 



In 1828 Robert Brown reported his observation on the rapid and irregular motion of 
pollens suspended in water [ij. The phenomenon was explained theoretically later by the 
combinational effects of a fluctuating force and a viscous drag force resulting from colli- 
sions with solvent molecules P|. Both the phenomena and the theory are immeasurably 
important in modern physics and continual development has been made for centuries 3|. 
With the advent of recently developed nanotechnology and inspiration from the biological 
world, people now are interested in constructing small devices working at molecular level 
and activated by thermal noise. One example is the so-called Brownian ratchet Q], where 
particles drift unidirectionally over a switchable asymmetric and periodic potential profile. 
In this Letter we present multiwalled carbon nanotubes based nanomachines, which can be 
fuled by thermal fluctuation. With appropriate control, it can act as one-dimensional device 
working at a length and time scale of nano. 

Carbon nanotubes, composed by one or more layers of cylindrical graphene shells, at- 
tracted much interest from various disciplines recently . These intriguing quasi-one dimen- 
sional structures possess outstanding material properties, such as high aspect ratio, chemical 
and thermal stability, high modulus and strength, and good transport properties. It also 
holds the promise for the building blocks in the next-generation nano-electromechanical 
systems (NEMs) [6]. Except for their unique mechanical and electronic properties, one of 
the most fascinating features of this novel material is that within the graphene shells, the 
sp^ bonds between carbon atoms are the strongest bond in nature, while between different 
shells only very weak van der Waals interactions. This high anisotropy provides the carbon 
nanotubes based devices with great intershell mobilities without losing structure stability 
as operating. 

Especially, model devices, such as oscillators and actuators based on multiwalled carbon 
nanotubes, were proposed Q and experimentally implemented in past few years. These 
devices, where graphene shells act as moving parts, have many impressive properties such 
as ultra high frequency and atomic precision. In the axial oscillator model 3, 1^ , the inner 
core within a open-end holder was pulled out at first and released then. The intershell 
restoring van der Waals force would drag the core back into the holder. Because of its 
ine..t,a a.d the .eg.,g,b.e .es,sta,.e fo.ce fa. the co.e w.ll sUde bae. a„d fo.K afte™a,.ds. 
The operating frequency was predicted to be up to gigahertz (GHz). Because of its great 
potential in engineering applications, lots of studies have been carried out on the dynamics. 
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mainly through atomistic simulations These works confirmed the high 

frequency on the order of GHz, but the energy dissipation was found to be rather severely. 
The axial oscillation initialized by extrusion decayed quickly in few nanoseconds. It's far 
away from our intuition and can not be explained by the previously predicted small static 
friction force. 

In fact, in such a device working at high speed on the order of 1.0 nm/ps, sliding will dissi- 
Eateb^^aMng kinetic eue.gy to othe. i„te..a. .nodes such as bending and .adial breathing 



15|,ll6|. The rate of leaking was found to depend on the temperature, chiralities and the 



system size. More importantly, open tube ends in the finite systems also make significant 
contribution jisl- Thus from an application point of view, low dissipation can be expected by 
lowering commensurabilities and avoiding the ends crossing of nanotubes during operating. 

To verify this, we study the incommensurate double-walled carbon nanotubes (DWNTs) 
(7,7)@(21,^as shown in Figure 1(a) by molecular dynamics simulation jl8]. The Dreiding 
forcefield |l2| is used to describe the interatomic interaction. Firstly the intertube van der 
Waals potential energy profile as a function of the axial displacement is illustrated in Figure 
1(a). It can be seen that the core is confined within a steep well of V^nds = 6.47 meV per 
atom provided by the tube ends repulsion. However inside the well, the local corrugation 
against sliding is much smaller, about V^oc = 0.009 meV per atom, which is comparable 
to the thermal energy at room temperature. So the dissipation of sliding inside the ends 
potential energy wells should be small by avoiding the ends effect following the arguments on 
the dissipation mechanism given above. In our simulation the system is firstly coupled with 
a 100 K bath for 200 ps. After the thermal equilibria, an axial velocity is assigned to the 
inner tube and the nanotubes start to oscillate between the outer tube ends. The relative 
displacement evolution of the DWNTs is shown in Figure 1(b). The axial oscillation is found 



to have negligib 
those reported 



e dissipation during the 5 nanosecond simulation. This result contradicts 
in which the length of the DWNTs are comparable and the 



inner tube is pulled out of the holder, so that tube ends cross and instabilities are involved 
The ultra-low dissipation observed here provides possibilities for energy supplement 
and experimental observation with a much larger timescale than nanoseconds. 

Physically speaking, the dissipation of sliding is the redistributing process from the sliding 
to internal modes. When the system reaches equilibrium, each mode, including axial sliding, 
holds a portion of the kinetic energy ^ksT following the equipartition theorem. And the 
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FIG. 1: (Color online) (a) Atomic structure and axial potential energy profile of the DWNTs 
(7,7)@(21,0), with a typical intershell distance of 0.34 nm. The inner core (red), with a length of 
2.5 nm, is confined by the ends repulsion potential well Vends = 6.47 meV per atom. While inside 
the well, there are three orders smaller corrugation T^oc = 0.009 meV per atom provided by the 
interaction with the 20 nm long outer tube (gray), as shown in the inset, (b) Axial sliding of the 
DWNTs, initialized by an axial speed at lOOK. The oscillation amplitude does not decay during 
the 5 nanoseconds simulation, which suggests ultra-low dissipation of the system. 

DWNTs will oscillate, in a fluctuative manner, relative to each other. This kind of fluctuative 
o.c,..a«o„ was observed by Se™,.ie a.d Gaspard .ecenUy Q. Howeve. a. .e„ti„„ed 

above, because the ends effects are involved in their work, the local potential wells at the 
equilibrium state are too steep, so the fluctuative the position (r^)^ = ^Jk^T /k is bounded 
in few Angstroms, where k is the effective elastic constant of the well. In contrast, the axial 
barrier in our model is negligibly small here, so large amplitude diffusion of the inner core 
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could be expected. Considering the intershell sliding motion as a normal mode, and the 
reduced mass as u = , where mi„ and mout are the mass of the inner and outer 

tube respectively. The velocity needed to overcome the local barrier V^oc is estimated to be 
vioc ~ 0.05 nm/ps. As the axial speed excited by the thermal noise exceeds this critical 
value, diffusion over the local barriers will be established. 

For this purpose, the DWNTs are investigated under thermal coupling. In the simulations, 
the outer tube is connected to Nose-Hoover thermostats at various temperature from 20 K 
to 300 K, and the core feels the thermal coupling through collision with the outer tube 
only. After equilibrating time of 200 ps, the data is acquired as integrating the equations of 
motion and plotted in Figure 2. At lower temperature, 20 K for example, the sliding motion 
is trapped within the 0.2 nm wide local barriers V^oc- The axial speed fluctuates with a 
peak value of 0.01 nm/ps which is much less than fioc. However at higher temperature 
as 100 K and 300 K, where more kinetic energy distributed on the sliding mode, global 
diffusion and oscillation between the tube ends is excited. One remarkable feature of this 
thermally excited oscillation is its large amplitude extending to the whole outer tube and 
high frequency on the order of GHz, which can be further tuned by the design of DWNTs. 

The oscillation behavior observed can be resolved into two regimes: (1) the bounded 
fluctuation state (BFS) (2) large-scale diffusion extending to the whole outer tube. In 
the BFS, axial speed fluctuates but is lower than the critical value fioc- Specifically, for 
the thermally excited oscillation at 100 K, kinetic energy distributed on sliding from the 
equipartition theorem is ^k^T, corresponding to average axial speed of 0.016 nm/ps, which is 
consistent with the simulation results. This consistency suggests the equilibrium fluctuation 
nature of the BFS. When the axial speed exceeds f loc at some moment, the inner tube jumps 
out of the local barriers and diffuse over them. Because of the ultra-low sliding dissipation 
of the DWNTs, the diffusion speed almost keeps constant at a value close to fioc. 

In addition, arguments should be addressed that the sliding fueled by the thermal bath 
observed in the DWNTs is not a violation of the second law in thermodynamics. The 
oscillational behavior is fluctuative and the axial speeds, the trapped time and the thus 
the oscillation frequency are not deterministic. For application, external controls such as 
external field or temperature gradient which introducing energy, must be applied to 
obtain a regular and controllable motion. 

Besides of the axial diffusion, angular motion of the tubes is also observed to be excited 
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FIG. 2: (Color online) (a) Axial displacement of the DWNTs coupling with a thermal bath at 20 K 
(red), 100 K (green), and 300 K (blue) respectively. The fluctuative oscillation at low temperature 
is confined in the local potential barrier of 0.2nm as shown in Figure 1 (a), while in temperature 
higher than 100 K large amplitude axial sliding is excited. Besides the constant speed sliding over 
local barriers, there are also some flat regions in the displacement curve. In these regions the 
axial diffusion is temporarily trapped, resulting from the thermally induced structure distortion as 
explained in the text later, (b) Corresponding axial velocity profile. There exists a critical value 
of the sliding speed Vct, values higher than that correspond to the linear diffusion in (a), and lower 
values correspond to the trapped region. The critical value Vcr = 0.05 nm/ps, read from the curve, 
coinciding well to the estimated velocity vioc to overcome the local potential barrier. 
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FIG. 3: Relative angular velocity of the double- walled nanotube, coupled with thermal bath at 
100 K. The velocity fluctuates with a mean value of 0.04 rad/ps. 

in the simulations. In contrast to the axial oscillation, unidirectional rotation is driven with 
average speed at 0.04 rad/ps. However, the amplitude is also fluctuative as shown in Figure 
3. 

Finally we consider the irregular parts of the oscillation profile shown in Figure 2(a), 
especially the flat regions in the curve of sliding amplitude at 100 K. The atomic structure 
of the nanotubes at the corresponding simulation time is extracted and shown in Figure 4. 
Structure distortion, which causes the outer tube to deviate from its cylindrical geometry, 
has been found. The local potential profile against sliding is calculated and plotted. It is 
found that the tube deformation induces an enhanced local barrier of 0.025 meV per atom, 
which requires a larger = 0.083 nm/ps. This value can not be exceed by the fluctuative 
speed shown in Figure 2. Thus the linear diffusion of the inner tube is trapped again and 
will not continue to slide until the local distortion is released. The distortion of graphene 
lattice results from the equilibrium fluctuation of other internal modes such as the bending 

n 

and rocking modes jl^ . Using multiwalled nanotubes or apply constraints on the outer tube 
may help to inhibit the deformation and avoid these instabilities. 

We have also performed simulations under higher temperature 1000 K and 1500 K. The 
distortion of the graphene shells are much more significant, and the global axial motion is 
frequently inhibited within few nanometers, although higher axial speeds over 0.1 nm/ps are 
excited. 
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FIG. 4: (Color online) The trapped state at 100 K corresponding to the flat region in displacement 
curve, at t = 4.8 ns, plotted in Figure 2(a). The cylindrical lattice of outer tube is found to be 
distorted. This imperfection causes an enhanced local well of 0.025 meV per atom against sliding, 
which requires a larger axial speed ( 0.083 nm/ps) to overcome. 

In conclusion, DWNTs based nanomachines driven by thermal bath are proposed in this 
Letter. Coupling with thermal bath at certain temperature, the inner tube can oscillate 
within the ends repulsion potential wells. The oscillation is a fluctuative phenomenon, with 
amplitude and operating frequency decided by the length of the outer tube. With external 
control applied, the nanomachine can be used to fabricate artificial nano devices. The physics 
behind this novel phenomenon is the kinetic energy hold by the translation or rotation 
at the thermal equilibrium is sufficient to overcome the corresponding potential barrier 
provided by the atomistically smooth potential surface of the graphene shell. Furthermore 
the ends potential provides a working space within which the severe dissipation by ends 
crossing is avoided. We have also investigated the confined dynamics of Ceo inside carbon 
nanotubes. Similar fiuctuative axial and angular motion are excited, furthermore the the 
potential profile can be tuned by varying the diameter and chirality of nanotube. Thus 
our observation suggests also the possibilities of fast molecule transport channel assisted by 
thermal noise and tuned by the tube geometry. On the other hand, this model system may 
also be theoretically interesting for studies of the nanoscale thermodynamics j^l • However, 
there still remain great but exciting experimental challenges for observing and controlling 
the nanoscale dynamics. 
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